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In this Chapter QCD interactions between a quark and an anti-quark 
are discussed. In the heavy quark limit these potentials can be related to 
quarkonia and l/m corrections can be systematically determined. Ex- 
citations of the ground state potential provide an entry point into the 
phenomenology of quark-gluon hybrids. The short-distance behaviour of 
non-perturbative potentials can serve as a test of resummation and con- 
vergence of perturbative expansions. Torelons and potentials between 
non-fundamental colour charges offer a window into the origin of the 
confinement mechanism and relate to effective string descriptions of low 
energy aspects of QCD. 



1.1 Motivation 

In order to avoid excessive overlap with the many introductions that al- 
ready exist in different places of this Volume, we shall elaborate on a very 
subjective motivation of studying interquark potentials in QCD, centred 
around the general theme of building bridges between models and QCD. 

QCD contains two related, non-perturbative features: the breaking of 
(approximate) chiral symmetry and the (effective) confinement of coloured 
objects such as quarks and gluons. Similarly, models of the QCD vacuum 
can roughly be divided into two classes: those that are primarily based 
on chiral symmetry and those that have confinement, for instance in the 
form of a confining potential, as their starting point. Often it is difficult 
to microscopically relate a particular model to QCD. For example, should 
the instantons that are evidenced in lattice simulations within a given pre- 
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2 Two Quark Potentials 

scription at a finite cut-off, be the same that are supposed to appear within 
instanton liquid models? The answer to this question is not known. 

There exist, however, two ways of systematically building bridges. Per- 
haps the most obvious one is the comparison of global properties like form 
factors, charge distributions in position space, potential energies or particle 
masses. Such a comparison circumvents the problem of identifying a one- 
to-one mapping of the degrees of freedom within a given model, which often 
might not qualify as a quantum field theory, onto objects that appear in the 
QCD vacuum. The model would then very much resemble the pragmatic 
use of analogies in the popular science literature: it is not that important 
to get things 100 % right if we cannot understand them 100 % anyway. 

A lot of purpose-engineered QCD information that is not directly acces- 
sible to experiment can be manufactured in lattice simulations. Moreover, 
in experiment there is only one world. On the lattice one has the freedom 
to vary the number of sea quark flavours n/, the quark masses, the number 
of gauge group colours, Nc, the volume etc. etc., away from the values that 
happen to be realized in nature for one reason or another. Sometimes (like 
in the case of quark masses) this is at present a necessity, in other cases it is 
pure virtue. In having not only one physical world but many virtual worlds 
at ones disposal, the applicability range and precision to be expected from 
any model can, in principle, be tested very stringently. Unfortunately, in 
practice, such interaction between lattice practitioners and model builders 
is still rather under-developed for various psychological, communication- 
related and dogmatic reasons and, in some cases, even out of fear, mistrust 
or over-confidence, one might speculate. 

The other way of building bridges is when a separation of scales occurs, 
in which the symmetries of QCD constrain the number of possible terms 
and allow for the systematic construction of an effective field theory. The 
most prominent examples are the chiral effective field theory (%EFT) that 
governs the low energy interactions of QCD as well as heavy quark effec- 
tive theory (HQET) and non-relativistic QCD (NRQCD) for heavy quark 
physics. While in %EFT the scale separation occurs through the sponta- 
neous breaking of the chiral symmetry by some collective gluonic effect, in 
HQET and NRQCD the scale separation is provided by the heavy quark 
mass to. Ideally, one would calculate the "high" energy Wilson coefficients 
of xEFT, such as the pion decay constant F n , in lattice simulations (or de- 
termine them from fits to experimental data). The low energy expansion, 
as a function of m ff , can be determined analytically (chiral perturbation 
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The Static QCD Potential 3 

theory). On the other hand, in HQET and NRQCD the low energy matrix 
elements can be provided by non-perturbative lattice simulations while the 
Wilson coefficients are calculable in perturbation theory, as long as m > A, 
where A denotes a typical non-perturbative scale of order 400 MeV. 

It is this latter heavy quark limit, in which the static QCD potentials 
that are discussed here — as well as in Chapter 4 — can be related to 
mesonic bound states and, in the case of three body potentials, baryonic 
bound states. In doing this, QCD can be reduced to non-relativistic quan- 
tum mechanics within this particular sector, with the help of reasonable 
assumptions which themselves can be tested in lattice simulations. For ex- 
ample, in the case of quarkonia, QCD itself tells us what "potential model" 
we have to choose. 

It is also possible to fit the spectra of light mesons and baryons, as- 
suming phenomenological potentials. These potentials cannot be related in 
a systematic way either to QCD or to static potentials as calculated from 
Wilson loops. However, in spite of some rather dubious assumptions that 
are implicitly folded into such non-relativistic or relativized quark-potential 
models, it can still be instructive to compare the parametrizations that are 
commonly used in this context with those that are relevant in quarkonium 
physics. After all, most physics appears to interpolate smoothly between 
hadrons made out of light and heavy quarks. There also appears to be 
an intimate connection between, on one hand, the broken chiral symmetry 
that is most relevant in the light hadronic sector — but plays little direct 
role for quarkonia with masses much larger than the chiral symmetry break- 
ing scale — and, on the other hand, the confining potential and flux tubes 
that seem to be the vacuum excitations that are responsible for interactions 
between the slowly moving heavy quark degrees of freedom. 



1.2 The Static QCD Potential 

The very definition of a "potential" requires the concept of "instantaneous" 
interactions: a test particle has to interact with the field induced by a source 
on a time scale short enough to guarantee that the relative distance remains 
unaffected. Whenever the relative speed of the two particles becomes rel- 
ativistic, the underlying assumption of a constant time difference between 
cause and effect is obviously violated: only in non-relativistic systems, i.e. 
as long as the typical interaction energies (E) within bound states are small 
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4 Two Quark Potentials 

compared to the particle masses (M), can we define a potential as a func- 
tion of coordinates such as the distance r, the spin S, angular momentum 
L and relative momentum p. While for interactions between elementary 
charges in QED this is always the case as E/M = O(ap) <C 1, QCD im- 
plies typical binding energies of order 400 MeV and only the bottom, and 
eventually the charm quark, can be regarded as non-relativistic. 

Another example in which the non-relativistic approximation is justified 
are nucleon-nucleon potentials, Vnn (0 ■ Although there exist attempts to 
extract this information also from QCD, by employing lattice simulations 
(cf. Chapter 4 Subsec. 4.4.2), it is fair to predict that quantitatively re- 
liable information, using the presently available methodology, might not 
become available within this decade. However, in this case a wealth of phe- 
nomenological information exists from experiment. For the quark- antiquark 
potential it is exactly the opposite. All experimental information is model 
dependent and rather indirect as quarks never appear as free particles in 
nature. However, this potential is among the most precisely determined 
quantities that have been calculated so far on the lattice 

As a starting point one can make the test charges infinitely heavy, pro- 
hibiting any change in their relative speed and study the "static" limit. 
To this end we shall introduce the Wegner- Wilson loop and derive its re- 
lation to the static potential. Subsequently, expectations of this potential 
from exact considerations, strong coupling and string arguments as well as 
from perturbation theory are presented. Lattice results are then reviewed. 
Finally, the model is extended to non-static quarks and the form of the 
resultant potential in coordinate space is compared with its counterpart 
(w-meson exchange) in the nucleon-nucleon interaction. 



1.2.1 Wilson loops 

We will derive the relationship between the expectation values of Wegner- 
Wilson loops and the potential energy V(r) = Eo(r) between two colour 
charges, separated by a distance r. This is a technical but instructive 
exercise. The final result is displayed in Eq. 11.171 

The Wegner- Wilson loop was originally introduced by Wegner |T] as an 
order parameter in Zi gauge theory. It is defined as the trace of the prod- 
uct of gauge variables U x ^ along a closed oriented contour 5C, enclosing 
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an area C, 




While the loop, determined on a gauge configuration* {U x p}, is in gen- 
eral complex, its expectation value is real, due to charge invariance: in 
Euclidean space we have (W(C)) = (W*(C)) = (W(C)>*. It is straight 
forward to generalise the above Wilson loop to any non-fundamental repre- 
sentation D of the gauge field, just by replacing the variables U x ,y, with the 
corresponding links U®^. The arguments below, relating the Wilson loop 
to the potential energy of static sources, go through independent of the 
representation according to which the sources transform under local gauge 
transformations. In what follows, we will denote a Wilson loop, enclosing a 
rectangular contour with one purely spatial distance, r, and one temporal 
separation, t, by W(r,t). Examples of Wilson loops on a lattice for two 
different choices of contours 6C are displayed in Figure 11.11 




Fig. 1.1 Examples of rectangular on- and off-axis Wilson loops with temporal extent 
t = 5a and spatial extents n = 3a and T2 = 2\/2a, respectively. 

*U \i <E SU(3) a gauge group element, pointing into direction fi £ {1,2,3,4}, located 
between positions x and x + a/t, where a is the lattice spacing and i is a lattice point 
within the (finite) 4-volume. 
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6 Two Quark Potentials 

In Wilson's original work [2], the Wilson loop had been related to the 
potential energy of a pair of static colour sources by using transfer ma- 
trix arguments. However, it took a few years until Brown and Weisberger 
attempted to derive the connection between the Wilson loop and the ef- 
fective potential between heavy, not necessarily static, quarks in a mesonic 
bound state 3 . Later on, mass dependent corrections to the static poten- 
tial have been derived along similar lines [41 5| and the approach has been 
made systematic within the framework of EFTs [61 71819) . In Sec. 11.31 we 
will discuss these developments in more detail. Here, we derive the connec- 
tion between a Wilson loop and the static potential between colour sources 
which highlights similarities with the situation in classical electrodynamics. 

For this purpose we start from the Euclidean Yang-Mills action, 



S =1^2 l^xF^. (1.2) 



The canonically conjugated momentum to the field, Af, is given by the 
functional derivative, 

5{diA°;) g 2 g 

The anti-symmetry of the field strength tensor implies — 0. In order to 
obtain a Hamiltonian formulation of the gauge theory, we fix the temporal 
gauge i.e. A\ = 0. In infinite volume such gauges can always be found. On 
a toroidal lattice this is possible up to one time slice t' , which we demand 
to be outside of the Wilson loop contour, i.e. t' > t. 
The canonically conjugated momentum, 

'6 A" ' 



< = -i^-r-, (1-4) 

/i £ An. ' V ' 



now fulfils the usual commutation relations, 

L4>£] = z^ b , (1.5) 
and we can construct the Hamiltonian, 

H = Jd 3 x (rfdtAI - ^F^F^j =\J d 3 x (EfE° - BfBf) (1.6) 

that acts onto states ^[^4^]. In the Euclidean metric the magnetic contribu- 
tion to the total energy is negative. Note that we can also add a fermionic 
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term Y^,f Q.f\rf^n + m /]9/ to the action above. In this case the momentum 
conjugate to the Dirac spinor field is given by —iS/SqJ = 9/74, where 
a = 1, • ■ • ,Nc runs over the colour in the fundamental representation. 
Here / denotes the quark flavour and to/ the respective mass. 

A gauge transformation f2 can be represented as a bundle of SU(Nc) 
matrices in some representation R, CIr(x) = e luJ ( x ) t h. We wish to de- 
rive the operator representation of the group generators Tjj, that acts on 
the Hilbert space of wave functionals. For this purpose we start from the 
definition 



l + i 



|d 3 x W Q (x)T^(x) + 



* = * + (S*. (1.7) 



One easily verifies that SAt = Af — At = — (cW + i[Aj,cj]) = Djw(x). Wc 
then obtain 

= -- fA w a (x)(A£,r(x)^ (1.8) 
9 J 

where we have performed a partial integration and have made use of the 
equivalence 

Si , 
JT =--E i 1.9 
oAi g 

of Eos. II. 31 and II. 41 Hence we obtain the representation 

T a R = —{DiEiY, (1.10) 

i.e. the covariant divergence of the electric field operator is the generator 
of gauge transformations! Again note that had we included sea quarks into 
the action, we would have encountered an additional term — - Y]f Qfj4T a qf 
on the right hand side of this equation, where the generator is to be taken 
in the fundamental representation. It is trivial to generalize the equations 
below accordingly and the physical meaning is clear too: the vacuum has 
an intrinsic charge density distribution due to sea quarks. This then in turn 
allows for string breaking of the static potential. 
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8 Two Quark Potentials 

Let us assume that the wave functional is a singlet under gauge trans- 
formations = ^[Afj]. This implies that 

(A^r* = o, (l.n) 

which is Gauss' law in the absence of sources: ^ lies in the eigenspace of 
DiEi that corresponds to the eigenvalue zero. Let us next place an external 
source in the fundamental representation of the colour group at position r. 
In this case, the associated wave functional $ a , a = 1, • • ■ , Nc transforms 
in a non-trivial way, namely 

[R(Sl)*] a = Otftf/J. (1.12) 

This implies that 

(A£i)°* = -gS 3 (r)T a V, (1.13) 

which again resembles Gauss' law, this time for a point-like colour charge 
at position^ r. For non-fundamental representations D, Eq. 11.131 remains 
valid under the replacement T a — > Tg. 

Let us now place a fundamental source at position and an anti-source 
at position r. The wave functional \l/ r , which is an Nc X Nc matrix in 
colour space will transform according to 

= n aT (0)O^(r)* r , 75 . (1.14) 

One object with the correct transformation property is a gauge transporter 
(Schwinger line) from to r, 



*r = ^=C/ f (r,t) ' 



fNc VNc 



exp J dxA(x, t) 



(1.15) 



which on the lattice corresponds to the ordered product of link variables 
along a spatial connection between the two points. Since we are in the 
temporal gauge, A${x) — 0, the correlation function between two such lines 
at time-like separation t is the Wilson loop 

(W(r,t)) = ^<EMr,t)L$ a (r,0)>, (1.16) 

tOf course, on a torus, such a state cannot be constructed. Note that in our Euclidean 
space-time conventions Gauss' law reads [DE]°(x) = — p a (x), where p denotes the 
charge density. Again note that, in general, p a will automatically contain a contribution 
gqff4T a qf for each sea quark flavour /. 
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which, being a gauge invariant object, will give the same result in any gauge. 
Other choices of ^r, e.<7- linear combinations of spatial gauge transporters 
connecting with r, define generalised (or smeared) Wilson loops, Wy(r, t). 

We insert a complete set of transfer matrix eigenstates, |$ r .„), within 
the sector of the Hilbert space that corresponds to a charge and anti-charge 
in the fundamental representation at distance r, and expect the Wilson loop 
in the limit of large temporal lattice extent, L T a 3> t, to behave like 



where the normalisation convention is such that = (^xfr) — 1 and 

the completeness of eigenstates implies J2 n K^nl^lO)! 2 = 1- Note that no 
disconnected part has to be subtracted from the correlation function since 
^ r is distinguished from the vacuum state by its colour indices. E n (r) 
denote the energy levels. The ground state contribution E (r) — that will 
dominate in the limit of large t — can now be identified with the static 
potential T^(r), which we have been aiming to calculate. 

The gauge transformation properties of the colour state discussed above, 
which determine the colour group representation of the static sources and 
their separation r, do not yet completely determine the state in question: 
the sources will be connected by an elongated chromo-electric flux tube. 
This vortex can, for instance, be in a rotational state with spin A / 
about the inter-source axis. Moreover, under interchange of the ends the 
state can transform evenly (r) = g) or oddly (j) = u), where r\ denotes the 
combined CP parity. Finally, in the case of the one-dimensional A = 
representations, it can transform symmetrically or anti-symmetrically un- 
der reflections with respect to a plane containing the sources (a v — ±). It 
is possible to single out sectors within a given irreducible representation 
of the relevant cylindrical symmetry group [TJJ) Doo/i, with an adequate 
choice of ^f. A straight line connection between the sources corresponds to 
the Doah quantum numbers £+, where A = 0, 1, 2, • • • is replaced by capital 
Greek letters, S, II, A, • • • . Any static potential that is different from the 

ground state will be referred to as a "hybrid" potential. * Since these 
potentials are gluonic excitations they can be thought of as being hybrids 
between pure "glueballs" and a pure static-static state; indeed, high hybrid 
excitations are unstable and will decay into lower lying potentials via the 

tCompare the discussion in Sec. 2.3.1 of Chapter 2 




(1.17) 



n 
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10 Two Quark Potentials 

radiation of glueballs. 



1.2.2 Exact results 

We identify the static potential V(r) with the ground state energy Eq(t) 
of Eq. 11.171 that can be extracted from the Wilson loop of Eq. 11.11 By ex- 
ploiting the symmetry of a Wilson loop under an interchange of the space 
and time directions, it can be proven that the static potential cannot rise 
faster than linearly as a function of the distance r in the limit r — > oo jllj . 
Moreover, reflection positivity of Euclidean n-point functions |12l 13) im- 
plies convexity of the static potential |14| . i.e. 

V"(r) < 0. (1.18) 

The proof also applies to ground state potentials between sources in non- 
fundamental representations. However, it does not apply to hybrid excita- 
tions, since in this case the required creation operator extends into spatial 
directions, orthogonal to the direction of r. Due to positivity, the potential 
is bound from below. § Therefore, convexity implies that V(r) is a mono- 
tonically rising function of r, i.e. 

V'(r)>0. (1.19) 

In Ref. |15|. which in fact preceded Ref. [TJ, somewhat more strict 
upper and lower limits on Wilson loops, calculated on a lattice, have been 
derived: Let a a and a T be temporal and spatial lattice resolutions. The 
main result for rectangular Wilson loops in representation D and d space- 
time dimensions then is 

(W(a a , a T )) rt ^ a " a ^ < (W(r, t)) < (1 - c y/*r+t/°r-* t (1.20) 

with c = exp[— 4(d— \)Dj3\. The resulting bounds on V{r) for r > a a read 

-ln(l-c) < a T V(r) < -—\n(W{a a ,a T )); (1.21) 

a a 

§The potential that is determined from Wilson loops depends on the lattice cut-off, a, and 
can be factorised into a finite potential V(r) and a (positive) self energy contribution: 
V(r;a) = V(r) + V so if(a). The latter diverges in the continuum limit (see Sec. 11.2151 . 
whereas the potential V (r) will become negative at small distances. Thus V (r; a) is 
indeed bounded from below by V(0)=0. 
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in consistency with Ref. |11|. the potential (measured in lattice units a T ) 
is bound from above by a linear function of r and it takes positive values 
everywhere. 

1.2.3 Strong coupling expansions 

Expectation values can be approximated by expanding the exponential of 
the lattice action, in terms of the inverse coupling (3 = 2Nc/g 2 , giving 
exp(— (3S) = 1 — [3S + • • • . This strong coupling expansion is similar to 
a high temperature expansion in statistical mechanics. When the Wilson 
action is used each factor (3 is accompanied by a plaquette and certain 
diagrammatic rules can be derived |21 161 171 181 19| . Let us consider a strong 
coupling expansion of the Wilson loop, Eq. 11.11 Since the integral over a 
single group element vanishes, 

dUU = 0, (1.22) 

to zeroth order, we have (W) = 0. To the next order in j3, it becomes 
possible to cancel the link variables on the contour 6C of the Wilson loop 
by tiling the whole minimal enclosed (lattice) surface C with plaquettes. 
Hence one obtains the expectation value |18l 201 211 22| 

{W{C)) ~ { \fi/2N$\ - alCaiSC) + ■■■, N C >2} (1 ' 23) 

for SU(Nc) gauge theory. If we now consider the case of a rectangular 
Wilson loop that extends r/a lattice points into a spatial and t/a points 
into the temporal direction, we find the area law, 

(W(r,t))=exp[-a d rt]+-- - , (1.24) 

with a string tension 

a d a 2 = -d r In — . (1-25) 
lo 

The numerical value of the denominator applies to SU(3) gauge theory; 
the potential is linear with slope ad, and so colour sources are confined 
at strong coupling. Here d r — (\ri\ + 1 + I^DA > 1 denotes the ratio 
between lattice and continuum norms and deviates from d r = 1 for source 
separations r that are not parallel to a lattice axis. The string tension 
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of Eq. 11.251 depends on d r and, therefore, on the lattice direction; 0(3) 
rotational symmetry is explicitly broken down to the cubic subgroup Oh. 
The extent of violation will eventually be reduced as one increases (3 and 
considers higher orders of the expansion. Such high order strong coupling 
expansions have indeed been performed for Wilson loops and glueball 
masses |24| . Unlike standard perturbation theory, whose convergence is 
known to be at best asymptotic (251 2fi| . the strong coupling expansion is 
analytic around (3 = j2Z] and, therefore, has a finite radius of convergence. 

Strong coupling SU (3) gauge theory results seem to converge for (3 < 5 
— see, for example, Ref. [20] ■ One would have hoped to eventually identify 
a crossover region of finite extent between the validity regions of the strong 
and weak coupling expansions |28j . or at least a transition point between 
the leading order strong coupling behaviour, a 2 oc — ln(/3/18), of Eo. 11.251 
and the weak coupling limit, a 2 oc exp[— 27r/3/(3/3o)], set by the asymptotic 
freedom of the QCD /3-function. However, even after re-summing the strong 
coupling series in terms of improved expansion parameters and applying 
sophisticated Pade approximation techniques |29| . nowadays such a direct 
crossover region does not appear to exist, necessitating one to employ Monte 
Carlo simulation techniques. In Fig. 11.21 we compare the strong coupling 
expansion for the string tension, calculated to 0(g~ 2i ) i.e. 0{(3 12 ) 
with results from lattice simulations. While at large (3 the lattice results 
approach the weak coupling limit, there appears to be no overlap between 
weak and strong coupling and neither between strong coupling and lattice 
results in the region of interest {(3 ~ 6 corresponds to a lattice spacing 
a -1 « 2 GeV). We have taken the np = value of the QCD A parameter 
as determined non-perturbatively in Ref. |30| as normalization for the weak 
coupling expansion. The error band of the 0(g 6 ) expectation is due to the 
corresponding statistical uncertainty. The 0{g A ) central value lies within 
this band. There is no normalization ambiguity in the strong coupling 
results. Also at (3 > 5, the quality of convergence of the strong coupling 
expansion diminishes. This break-down might be related to a roughening 
transition that is discussed, for example, in Refs. |31I32| . 

We would like to remark that the area law of Ea. ll.24| is a rather general 
result for strong coupling expansions in the fundamental representation of 
compact gauge groups. In particular, it applies also to U(l) gauge theory 
which we do not expect to confine in the continuum. In fact, based on du- 
ality arguments, Banks, Myerson and Kogut |33j have succeeded in proving 
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the existence of a confining phase in the four-dimensional theory and sug- 
gested the existence of a phase transition while Guth [j^] has proven that, 
at least in the non-compact formulation of U(l), a Coulomb phase exists. 
Indeed, in numerical simulations of (compact) U (1) lattice gauge theory two 
such distinct phases were found [351 3f>| , a Coulomb phase at weak coupling 
and a confining phase at strong coupling. The question whether the con- 
finement one finds in SU(Nc) gauge theories in the strong coupling limit 
survives the continuum limit, p — > oo, can at present only be answered by 
means of numerical simulation (and has been answered positively). 




5 5.5 6 6.5 7 

P 

Fig. 1.2 Comparison of strong and weak coupling expansions with non-perturbative 
results for the string tension. The upper lines correspond to the strong coupling limit 
beginning with ln3cj 2 (dashed line) — the lowest order term from Eg. II. 2^)1 for 5(7(3). 
The other upper lines show the additional approximations up to 0(g~ 24 ). The lower 
lines correspond to the weak coupling limit up to 0(g a ) 
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1.2.4 String picture 

The infra-red properties of QCD might be reproduced by effective theories 
of interacting strings. String models share many aspects with the strong 
coupling expansion. Originally, the string picture of confinement had been 
discussed by Kogut and Susskind [SZj as the strong coupling limit of the 
Hamiltonian formulation of lattice QCD. The strong coupling expansion of a 
Wilson loop can be cast into a sum of weighted random deformations of the 
minimal area world sheet. This sum can then be interpreted to represent 
a vibrating string. The physical picture behind such an effective string 
description is that of the electric flux between two colour sources being 
squeezed into a thin, effectively one-dimensional, flux tube or Abrikosov- 
Nielsen-Olesen (ANO) vortex [381 391 41)1 41| . As a consequence, this yields 
a constant energy density per unit length and a static potential that is 
linearly rising as a function of the distance. 

One can study the spectrum of such a vibrating string in simple mod- 
els |31I42I 45] . Of course, the string action is not a priori known. The 
simplest possible assumption, employed in the above references, is that the 
string is described by the Nambu-Goto action (441 4T>] in terms of d — 2 
free bosonic fields associated to the transverse degrees of freedom of the 
string. In this picture, the static potential is (up to a constant term) given 

by mm 



(d-2W (d-2)n (d-2) 2 ir 2 

V(r) = ar\l-± ir- = ar-^ '- ^ . 1.26) 

1 ' V I2ar 2 24r 1152crr 3 V ' 

For a fermionic string 07] one would expect the coefficient of the correction 
term to the linear behaviour to be only one quarter as big as the Nambu- 
Goto one above. In the bosonic string picture, excited levels are separated 
from the ground state by 



V 2 (r) = V 2 (r) + (d- 2)ima = 



+ (d-2)nn 
2r 



(1.27) 



with n assuming integer values. It is clear from Eq. 061 that the string 
picture at best applies to distances 



For d = 4 one obtains r c « 0.33 fm, when using the value y/a w 430 MeV 
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from the p, 0,%, ■ ■ ■ Regge trajectory. 

The expectation of Eq. 11.261 has been very accurately reproduced in 
numerical simulations of Z-i gauge theory in d = 3 space-time dimensions 
[481 49| . In contrast, for d = 4 SU(3) gauge theory the spectrum of hybrid 
potentials still differs significantly from the expectation of Eq. 11.271 for 
distances as large as 2 fm [501 5 1| . However, qualitatively the string picture 
is supported by the £17(3) data too, since the hybrid potentials at large r 
are found to group themselves into various bands that are separated by 
approximately equidistant gaps. These will eventually converge to values 
7r/r, at even larger distances than accessible at present. Such an observation 
would support the existence of a bosonic string description of confining 
gauge theories in the very low energy regime |52l 531 541551 56) . Of course, 
in d < 26, the string Lagrangian is not renormalisable — but only effective 
— and higher order correction terms like torsion and rigidity will in general 
have to be added |57|. 

It is hard to disentangle in d — 4 the (large distance) 1 jr term, expected 
from string vibrations, from the perturbative Coulomb term at short dis- 
tances. However, a high precision attempt has been made recently, with 
promising results |58| . As an alternative, three-dimensional investigations 
(where perturbation theory yields a logarithmic contribution) have been 
suggested |59I60| . Another way out is to determine the mass of a closed 
string, encircling a boundary of the lattice with a spatial extent I = L a a 
(a torelon) , which is not polluted by a perturbative tail. The bosonic string 
expectation in this case would be |59| 

E n (l)=al- K ' +■•■■ 1.29 
6 / 

The naive range of validity of the picture is I ^> l c — 2r c « 0.66 fm. The 
numerical value applies to d = 4 from Eq. 11.281 An investigation of the 
finite size dependence of the torelon mass in d = 4 SU(2) gauge theory has 
been performed some time ago by Michael and Stephenson |61j who found 
excellent agreement — on the 3 % level — with the bosonic string picture 
already for distances lfm < I < 2.4 fm, quite close to l c . Qualitative 
agreement has also been reported by Teper 62j from simulations of SU(2), 
SU(3), SU(A) and SU(5) gauge theories in three dimensions as well as 
in a recent study of four-dimensional SU(2N) gauge theories by Lucini et 
al. 

The bosonic string picture prediction of the free energy, calculated from 
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Polyakov line correlators, at finite temperatures T is similar to Eg. 11.291 



im<P*(r)P(0))=<r(/3)r + - 



CT ( T )= a Q T 



(1.30) 



with validity for r ^> T = aL T [HI]. The Polyakov line is defined as 



P(x) = Tr < T 



exp | z / dx4 Ai(x) 
o 



= Tr J] C/ X)4 



(1.31) 



where T denotes time ordering of the argument. The dependence of the 
effective string tension on the temperature has been checked for rather 
low T- 1 < 1.24 T" 1 « 0.93 fm in studies of SU{3) gauge theory |65I66| . 
Although the sign of the leading correction term to the zero temperature 
limit is correct, the difference comes out to be larger than predicted. It 
would be interesting to check whether the result will converge towards the 
string expectation at lower temperatures. 



1.2.5 The potential in perturbation theory 

The strong coupling expansion is specific to the lattice regularisation.^ 
However, the expectation value of a Wilson loop can also be approximated 
using standard perturbative techniques. 

We will discuss the leading order weak coupling result that corresponds 
to single gluon exchange between two static colour sources which, although 
we neglect the spin structure, we will call "quarks" for convenience. From 
the Lagrangian, Cym = ^TrF^i 7 ^, one can easily derive the propagator 
of a gluon with four-momentum q, 

GtM^a 2 ^: (1-32) 

where /z, v are Lorentz indices and a,b — 1, ■ ■ • Na label the colour gener- 
ators with N\ — Nq — 1 for SU(Nc)- The same calculation can be done 
starting from a lattice discretised action. The Wilson action yields the 
result of Eci. I1.32I ud to the replacement 

<ln ~+ = ^ sin (^r) • (° 3 ) 

^However, new strong coupling methods have been developed in the large Nq limit, 
based on the Maldacena conjecture of QFT/AdS correspondence. 
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The Static QGD Potential 17 

Other lattice actions yield slightly different results but they all approach 
Ea. 11.321 in the continuum limit, a — > 0. Momentum space potentials can 
be obtained from the on-shell static quark anti-quark scattering amplitude: 
the gluon interacts with two static external currents pointing into the pos- 
itive and negative time directions, A^ a/3 = S^T^ and A'^ s = -8 vA T^ s . 
Hence, we obtain the tree level interaction kernel 

K aPlS {q) = ~T^ S . (1.34) 

For sources in the fundamental representation, the Greek indices denote 
the colour indices of the external currents running from 1 to Nc and the 
quark anti-quark state can be decomposed into two irreducible representa- 
tions of SU(N C ), 

N C ®N^ = 1®N A . (1.35) 

We can now either start from a singlet or an octet" initial <f>/3 7 = QpQ* 
state, 

*%, = <W (1-36) 



1 



where the normalisation is such that & % a p& 3 a a = 8 V . A contraction with 
the group generators of Eq. 11.341 yields 

^T^T; S = C F ^ aS , (1.38) 

^ A W S = -7^$^, (1.39) 

where Cp = Na/(2Nc) is the quadratic Casimir charge of the fundamental 
representation. 

We end up with the potentials in momentum space, 



governing interactions between fundamental charges coupled to a singlet 
and to an octet, respectively: the force in the singlet channel is attractive 
while that in the octet channel is repulsive and smaller in size. 

II We call the state an "octet" state, having the group SU(3) in mind. 
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How arc these potentials related to the static position space inter-quark 
potential, defined non-perturbatively through the Wilson loop, 

V(r) = - lim —ln<W(r,t))7 (1.41) 
t^oo at 

The quark anti-quark state creation operator, ^ r , within the Wilson loop 
contains a gauge transporter and couples to the gluonic degrees of freedom. 
Thus, in general, it will have overlap with both, QQ* singlet and octet 
channels**. Since the singlet channel is energetically preferred, i.e. V s < V a , 
we might expect the static potential to correspond to the singlet potential. 
Up to order g 6 this is indeed the case: to lowest order, the Wilson loop — 
defined by the closed contour SC — is given by the Gaussian integral 

(W(r,t))=exp|-i Jd 4 xd 4 yJ«(x)G a Jl(x~y)Jt(y)y (1.42) 

where J" = ±T a if (x,fx) G SC and = elsewhere 1 " 1 '. Eg. 1 1.421 implies 
for t » r 



(W(r,t)) = exp (c F g 2 t J ' dt 1 [G(r,t') - G(0,t') 



(1.43) 



We have omitted gluon exchanges between the spatial closures of the Wilson 
loop from the above formula. Up to order g 6 (two loops), such contribu- 
tions result in terms whose exponents are proportional to r and rjt and, 
therefore, do not affect the potential of Eg. 11.411 The propagator G a ^ v {x\ 
the Fourier transform of G°^ v (q) in Eq. 11.321 contains the function 

f d A a e iqx f°° 11 
G ^=/(4^' L dXiG{x)= T,-r ^ 



After performing the ^integration, we obtain 

V(v, n) = -C F — + V seU {fj,), (1.45) 
r 

**Of course, for a quark and anti-quark being at different spatial positions, the singlet- 
octet classification should be consumed with caution in a non-perturbative context. 

"Note that this formula, which automatically accounts for multi-photon exchanges, is 
exact in non-compact QED (excluding fermion loops) to any order of perturbation 
theory. However, in theories containing more complicated vertices, like non-Abelian 
gauge theories or compact lattice U(l) gauge theory, correction terms have to be added 
at higher orders in g 2 . 
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where a s = g 2 /(47r). The piece 

W = C F9 *J^^=C F a s 2 -», (1.46) 

that linearly diverges with the ultra-violet cut-off /i, results from self-inter- 
actions of the static (infinitely heavy) sources. Comparing Eas. 11.401 and 
II .451 we indeed find 

V(q) = V s (q), (1.47) 

where 




^ r V(v), V{v) = V{v^)-V s M. (1.48) 



This self-energy "problem" is well known on the lattice and has also received 
attention in continuum QCD, in the context of renormalon ambiguities in 
quark mass definitions |67l 68) . 

Note that while V s corresponds to the static potential, the perturba- 
tion theory relevant for hybrid excitations of the ground state potential 
corresponds to V p)9"| . 

At order a A s a class of diagrams appears in a perturbative calculation 
of the Wilson loop that results in contributions to the static potential that 
diverge logarithmically with the interaction time |7(J) . In Ref. |71|. within 
the framework of effective field theories, this effect has been related to ultra- 
soft gluons due to which an extra scale, V Q — V s , is generated. Moreover, 
a systematic procedure has been suggested to isolate and subtract such 
terms to obtain finite singlet and octet interaction potentials between heavy 
quarks. 

The logarithmic divergence is related to the fact that Eq. 11.431 con- 
tains an integration over the interaction time. For large times and any 
fixed distance r, Wilson loops will decay exponentially with t. However, 
the tree level propagator in position space is proportional to (r 2 + t 2 )" 1 , 
i.e. asymptotically decays like t~ 2 only. We notice that the integral receives 
significant contributions from the region of large t as demonstrated by the 
finite t ^> r tree level result 

-MW(r,t)) = -Q^t- (arctan- - — 

(1.49) 



In 1 



+ (r + t)V aeU . 
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The tree level lattice potential can easily be obtained by replacing 
by <7 M from Eq. 11.331 and (in the case of finite lattice volumes) the integrals 
by discrete sums over lattice momenta, 



2-77 rii 
Qi = 7 , 



L 
~2 



zz + i ... 



The lattice potential reads 
V(r) 

where 

1 
r 



V se \{(a) - C F a s 



iq-R 



4tt 



(1.50) 



(1.51) 



(1.52) 



and V S c\f(a) = Cfo. s [1/0]. We have neglected the zero mode contribution 
that is suppressed by the inverse volume {aL a )~^ . In the continuum limit, 
[1/r] approaches 1/r up to quadratic lattice artefacts whose coefficients 
depend on the direction of r while V se \[(a) diverges like 



F S eif(a) = Cpa^ 1 x 3.1759115 • • • 



(1.53) 



The numerical value applies to the limit, L a — > oo. Note that under the sub- 
stitution fi w 1.5879557 it /a, Ea. ll.53l is identical to Eg. II. 461 One loop com- 
putations of on-axis lattice Wilson loops can be found in Refs. [721 731 74*j . 
while off-axis separations in QCD with and without sea quarks have been 
realized in Ref. [7S]. The tree level form, Eg. 11.511 is often employed to 
parameterise lattice artefacts — see Subsec. 5.9.3.1 of Chapter 5. 



1.3 Quark- Antiquark Potentials between Non-Static Quarks 

** So far in this chapter the role of the interquark potential has been to test 
various limits and models for QCD. For example, in Fig. II. 21 a comparison 
was made between the weak and strong coupling limits of QCD and also 
with the corresponding lattice results. Likewise, in Subsec. 11.2.41 string 
models were compared with the infra-red properties of QCD. This role 

**At this point the author "ran-out-of-steam" and so, because of the pressure of time, 
the editor felt that this chapter should be rapidly concluded — a task carried out by 
the editor himself with the semi-approval of the author. 
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is in stark contrast to that of the Nucleon-Nucleon potential [V(NN)] in 
nuclear physics, where V(NN) is mainly used as a stepping stone to the 
understanding of multi-nucleon systems. To this end, the many parameters 
needed to define V(NN) are first adjusted, more or less freely, to fit two 
nucleon experimental data, with — in some cases — values being imposed 
from meson- nucleon data or theories. Unlike V(QQ), the NN-potential 
is not able to predict reliable quantitative information about the input 
parameters to its theory — the only exception perhaps being an estimate of 
the ttNN coupling constant [7E|- The reason why the form of V(NN) is so 
complicated compared with that so far discussed for V(QQ), i.e. essentially 
VqQ{r) = —e/r + cr, is because the latter is the interaction between two 
static quarks. If, however, we go away from this limit, then immediately the 
spins of the quarks begin to play a role resulting in forms similar to those 
encountered in V(NN). These potentials should be reliable for describing 
bb, be and cc states, since the b and c quarks are still sufficiently heavy 
4.5 and 1.2 GeV respectively) to not require a relativistic treatment. 
However, as soon as s, d, u quarks are involved, relativistic effects become 
important and even the whole concept of an interquark potential should be 
questioned. Furthermore,, the hope that this potential between two quarks 
can account for multiquark systems has yet to be justified. Even so, this 
has not deterred its use as an effective interaction — a topic discussed in 
Chapter 5. 

1.3.1 Radial form of V(QQ) 

What is the form of V{QQ)1 To this question there is no unique answer, 
since — as with meson-exchange models of V(NN) — forms depend on 
the theory from which the potential is derived and in which it should be 
utilized. A good example of this ambiguity is the momentum dependence 
of the potential. Even though the correct relativistic scattering equation 
for two quarks is the Bethe-Salpeter equation, for practical reasons, this 
needs to be simplified to, say, the semi-relativistic Blankenbecler-Sugar or 
non-relativistic Schrodinger equations — see Subsec. 5.1.2.2 of Chapter 5 
for a very brief discussion of this. If the basic Blankenbecler-Sugar equation 
is used directly then the appropriate potential contains relativistic factors 
of the form Eq = \J M 2 + p 2 , where p is the momentum of the quark with 
mass M. However, it is often convenient to expand these momentum factors 
in powers of p/M resulting in a potential appropriate for a Lippmann- 
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Schwinger or Schrodinger approach. 

Just as the interaction between two static quarks contained two distinct 
parts — i.e. the gluon exchange term Vg = — -^ct s /r and the confining term 
Vc = cr — so can the interaction between a heavy quark of mass M and 
an anti-quark of mass m be expressed as 

V(Mm) = Vg (Mm) + V c (Mm), (1.54) 

where 

V B( Mm )( -\ a ,r = I - ^(r), M .. m - ^S 12 



1 \M 2 + m 2 2 1 1 M 2 -m 2 

"2H» [ 2M 2 m 2 + Mm~\ L - S+ 8r^ M 2 m 2 -<T m ).L + - ■ ■ , (1.55) 

obtained by a non-relativistic reduction of the one-gluon-exchange mecha- 
nism and 

i -m * \ cM 2 + m 2 r „ cM 2 -m 2 , Nr . . 

V c (Mm) = cr-- T L.S+- (<7 M -a m ).L + - ■ • . 1.56 

r 4^4^771^ r 8M z m^ 

These expressions for pg(Mm) and Vc(Mm) do not include effects from 
expanding the E = y/ M 2 + p 2 terms mentioned above — see, for example, 
Ref. |8()j . Also here, for simplicity, they do not include spin- independent 
terms proportional to p 2 and <r 3 ' (r) , since their form becomes ambiguous 
when going from a momentum space to a coordinate space representation. 
It should be added that it is convenient to have VoiMm) in coordinate 
space, since it is hard to deal with the linearly rising term cr from Vc(Mm) 
in momentum space. In addition to this ambiguity there are three others 
that enter in practice: 

(1) A crucial term in Ea ll.55l is the hyperfine interaction proportional 
to (TM-Cm, which splits the energies of pseudoscalar and vector 
mesons. However, it is seen that it is proportional to S(r) — a 
radial form that needs to be regulated before it can be used in a 
wave-equation. This can be accomplished in a variety of ways , of 
which the most simple is to make the replacement 

(5 ( r )^47J ex P(- aV ) ( L5? ) 



Februa ry 4, 2008 14:16 WorldScientific/ws-b8-5x6-0 chap terlrv 



Quark- Antiquark Potentials between Non-Static Quarks 23 

as was done in the original work of Godfrey and Isgur [77]. Another 
approach is to consider this term as an effective interaction that is 
added in first order perturbation theory using those wavefunctions 
generated by the rest of V(QQ). In this case the overall constant is 
considered to be essentially a free parameter although reasonable 
estimates can be made from models involving instantons |78| . 

(2) In Ea. 11.451 a s is the quark-gluon running coupling "constant" i.e. 
it depends on momentum k. This can be parameterised in a variety 
of ways e.g. from Ref. jZH] 

where there are two parameters 

i) the gluon effective mass parameter m g 240 MeV and 

ii) the QCD scale parameter A ps 280 MeV. 

In Ref. [SUj this is combined with the regulator of <5(r) by the re- 
placement 

(1-59) 

where E M = ^M 2 + (fc 2 /4) and E m = ^m 2 + (fc 2 /4). 

(3) In the above it is implicitly assumed that Vc(Mm) is purely scalar. 
However, there are reasons — both theoretical and phenomenolog- 
ical — suggesting that there could be a sizeable vector component. 
For example, the spin-orbit splitting from Vg(Miti) is oc 1/r 3 , 
whereas that in Vc(Mm) is oc — 1/r. This indicates that the natural 
spin-orbit ordering of a coulomb-like potential should be inverted 
for high partial waves — a feature not seen experimentally or in 
lattice calculations. This is discussed in Subsec. 5.9.3.2 in Chap- 
ter 5. 

Lattice calculations are able to isolate several of the components in the 
potential V(Mm) in Eq. 11.541 For example, the leading confinement term 
Vc(Mm) = cr of Ea. ll.56| is shown in Fig. 11.31 Another point of interest in 
this figure is that there is no sign of the expected flattening at r/ro ~ 2.4, 
where it becomes energetically favourable to create two mesons — denoted 
by the horizontal band {51 . For other components of the potential the 
lattice data is less precise, but still there is reasonable agreement with the 
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Fig. 1.3 The linearly rising confinement as calculated on a lattice 1811 compared with 
V c (Mm) = cr in Eq. ITCT 



expectations from Eas. I1.55l and n~56l Probably the most extensive lattice 
study of the various forms contained in the potential appear in Ref. 

1.3.2 Comparison with the form of V(NN) 

For those more familiar with internucleon potentials it may be of interest 
to compare the forms that appear in V(NN) with those in V(QQ). This is, 
probably, best done by writing down explicitly the w-meson contribution 
to the familiar One-Boson-Exchange-Potential, since this meson is the one 
that resembles most closely a gluon both being vector particles independent 
of isospin/flavour. In V(NN), the w-meson contributes much of the short 
range repulsion and also has a strong Spin-orbit potential. The form given 
in Ref. [E2 is 



Vuj exp(— mr) 



1 



2M 2 



v2 exp(-mr) exp(-mr) ^ 2 



2M 2 



2 exp(— mr) 



(1 + ^1) 
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2.1 



m 

4M 2 



1 

fhr 



1 



(mr) 2 



exp(— mr) 3 1 1 rf 



exp(— mr) 



L.S, 



(1.60) 

where M is the nucleon mass and m — 783 MeV s» 4fm _1 is the mass of 
the to. When m —* to compare with one-gluon-exchange, this reduces to 



V(m 0) 



1 



2M 2 



V 2 - + -V 2 

r r 



1 



<r 1 .er 2 



1 s 3 



1 1 



i.5. 



(1.61) 



4M 2 r 3 1Z 2M 2 r 3 
This form is now quite similar to the one-gluon-exchange potential Vq in 
Ea. 11.551 — the only differences being the appearance of the spin- independent 
terms containing V 2 and S^(r), which are not unique in going from mo- 
mentum space to coordinate space. Of course, it is possible to stay through- 
out in momentum space, since — unlike Vc(Mm) — here there is no linearly 
rising potential cr. However, there are also differences in how Eqs. 11.551 
and 11.601 are treated. Normally in Eq. 11.601 the 5^ (r) terms do not ap- 
pear, since there are physically motivated models of form factors at the 
NNui- vertices, which smooth out such singularities. Also the results are 
not crucially dependent on these form factors since the main role of the 
w-meson is to generate a strong spin-isospin independent repulsion that 
essentially excludes the iViV-wavefunction from the region of the origin. 
In contrast, the one-gluon-exchange potential in Ea. 11.551 has an attractive 
spin-isospin independent term (— |a s /r). This means that uncertainties in 
the parametrizations in Eas. 11.571 and IT. 5 81 are more important. 



1.4 Conclusions 



This chapter has attempted to describe a few of the many topics that could 
be covered by the title. It has not aimed at being in anyway comprehensive 
— the selection being rather subjective with the following ideas in mind. 

1) The chapter began and ended with the theme that the interquark poten- 
tial V(QQ) could be viewed as a bridge for modelling QCD with somewhat 
the same role that the nucleon-nucleon potential V(NN) plays in nuclear 
physics. To further bring out the analogy between V(QQ) and V(NN) the 
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radial forms were compared and contrasted in Sec. 11.31 However, it should 
be added that this possible "interdisciplinary bridge" is not well understood 
and its study is rather neglected. 

2) The potential V{QQ) can be used for testing approximations to QCD 
and its lattice formulation. This is illustrated in Fig. 11.21 where it is shown 
that the weak and strong coupling limits do not overlap for the energy 
range of interest i.e. between 100 MeV and a few GeV. Furthermore, over 
this range neither limit agrees with the corresponding (non-perturbative) 
lattice calculation. This feature shows that we are "forced" into performing 
lattice calculations, since there seems to be no other way of treating QCD 
with the couplings of most interest. 

Here no mention has been made about deriving a discretized form of 
V(QQ). This topic has been the scene of much activity and comes under 
the heading of Non-Relativistic QCD (NRQCD). Here an expansion is made 
in terms of A/Mq with A ~ 1 GeV being a characteristic energy scale for 
non-perturbative effects. A few comments can be found in Subsec. 5.1.2.1 
of Chapter 5. 
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